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Abstract. We use the isotropy of the Cosmic Microwave Background to place 
stringent constraints on a possible electrical charge asymmetry of the universe. We 
find the excess charge per baryon to be q e - P < 10~ 26 e in the case of a uniform 
distribution of charge, where e is the charge of the electron. If the charge asymmetry 
is inhomogeneous, the constraints will depend on the spectral index, n, of the induced 
magnetic field and range from g e _ p < 5 x 10 -20 e (n = —2) to q e - p < 2 x 10 _26 e (n > 2). 
If one could further assume that the charge asymmetries of individual particle species 
are not anti-correlated so as to cancel, this would imply, for photons, q 7 < 10 _35 e; 
neutrinos, q v < 4x 10~ 35 e; and for heavy (light) dark matter particles <7dr 
(q dm < 4 x 10- 30 e). 
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1. Introduction 

With the substantial improvement of experimental and observational techniques in 
particle physics an astrophysics, it has become possible to test some of the assumptions 
that go into constructing models of fundamental interactions and the universe. The 
aim of this work is to place a cosmological constraint on the presence of an electrical 
charge asymmetry in the universe. The first cosmological analysis of a charged universe 
undertaken was in the context of the Steady State model of the universe by Bondi & 
Lyttleton in 1959 pQ, where an attempt was made to explain the recession of the galaxies 
through electromagnetic repulsion. To implement charge non-conservation, Maxwell's 
equations were modified to include a direct coupling to the vector potential, violating 
gauge invariance. Swann [2; confuted this analysis; Barnes jSj correctly reanalysed a 
similar model, using the Proca theory of electromagnetism. The spontaneous breaking of 
gauge symmetry, and the subsequent development of a charge imbalance, was considered 
by Ignatiev et al. [4 , and Dolgov & Silk p] used it as a mechanism of creation of 
primordial magnetic fields. An implementation of a homogeneous and isotropic charge 
density was proposed in the context of massive electrodynamics jH|, and the possibility 
of charge non-conservation has been analysed in brane world models [7], and in varying 
speed of light theories [H|. 

In principle, a potential charge asymmetry could be carried by a number of different, 
stable components in the universe. These can include: photons, neutrinos, dark matter, 
or a difference in the electron and proton charges. Experimental constraints, based on 
the lack of dispersion from pulsar signals, limit the photon charge to be g 7 < 5 x 10~ 30 e 
jU] (a less restrictive laboratory constraint, q 1 < 8 x 10~ 17 e, has also been established 
The analysis of the luminosity evolution of red giants in globular clusters leads 
to a constraint on the charge of neutrinos (or the charge of any sub-keV particle) of 
q v < 2 x 10 _14 e [H]. Direct, laboratory constraints based on gas-efflux measurements 
[12], electro-acoustic techniques [T3], and Millikan-type experiments using steel balls 
|14j . all limit the electron-proton charge difference to be q e - p < 10~ 21 e. For particles 
generally with masses less than ~ 1 MeV, Big Bang Nucleosynthesis rules out charges 
greater than 10~ 10 e jT^j. Concerning dark matter, constraints have been established on 
the fraction of dark matter due to charged heavy particles of the order of 10~ 5 [16J. 

In this work we will obtain general constraints on a possible charge asymmetry, 
characterised by a uniform or a stochastic distribution, placing the potential contribution 
of individual components in context. The premise will be that charge was generated at 
some very early time, but has been conserved in the period in which we are establishing 
constraints. 

We use Heaviside-Lorentz electromagnetic units with e = ^Aira, c = 1 = 87rG; 
Greek indices run from to 3 and Latin indices from 1 to 3. The scale factor is normalised 
to a(r ) = 1 today, where r denotes conformal time. 
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2. Model of a charged universe 

In this analysis we use the 1+3 covariant formalism developed in [17] : we consider a 
general class of homogeneous space-times for which it is possible to define a preferred 
velocity field, u a , determining the fundamental fluid flow lines and satisfying u a u a = — 1. 
The existence of this vector field generates a unique splitting of spacetime, given by the 
instantaneous three-dimensional rest space of an observer moving with 4-velocity u a 
and the one-dimensional space u a itself. The metric tensor of this spacetime, g a p can 
be written like g a p = h a/3 — u a up, where h a/3 is a projection tensor and represents 
effectively the spatial metric for the observer moving with u a . Any physical tensor 
field on this spacetime can be separated with respect to h a p and u a into space and 
time parts. The covariant derivative of a tensor field S af3 splits into a comoving time 
derivative S a/3 = u J S al3 ^ and a covariant spatial derivative D^S^ = h 1 s h a fl h f3 u S tJ-u -s- 
In particular, the splitting of the covariant derivative of the 4-velocity takes the form: 

u a -p = o- aP + u) af3 + ^Qh a p - u a u p , (1) 

where a a p is the shear tensor, u a p is the vorticity tensor, = 3d /a is the volume 
expansion, and u a = u a -,pw is the acceleration vector. The electromagnetic field tensor 
F a p = —Fp a is split into electric and magnetic fields as measured by an observer moving 
with u a : 

F a p = 2v,[ a Ep] + eap-yB" 7 , (2) 

where the electric and magnetic 3- vectors are given by E a = F a pvP B a = | e^i 7 ^ 7 , 
and e Q( g 7 = rj^su 5 is the volume element for the orthogonal rest space h a p. We define 
the current density 4- vector j a , where q = —j a u a is the charge density, and J a = h a p^ 
is the spatial current measured by u a . Then, Maxwell's equations F a ^-p = j a , 
F a /3;-y + Ffct;/3 + Fpj. a = become [TTJ [TH] 

D a E a + 2u a B a = q (3) 
D a B a - 2uj a E a = (4) 

KpE? - e a p y D?B^ = -J a + e^vPW + L a p + a a p - ^Qh a pj E 15 (5) 

hapB 13 + e a p 7 D^E^ = - e a p 7 u^EP + L a p + a a p - ^Qh a p^j B? , (6) 

where u a — h e^p^uj^ 1 is the vorticity vector. It appears that the motion of the observer 
affects the form of Maxwell's equations: beside the usual spatial divergences, curls and 
time derivatives of the fields, one has the appearance of extra terms representing the 
motion of the family of fundamental observers moving with 4-velocity u a . Moreover, 
the current conservation equation j a . a = takes the form 

q + Qq + D a J a + ii a J a = 0, (7) 

and the covariant form of Ohm's law is 

ja + U a Upj P = oF a pu p , (8) 
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where a denotes the conductivity. Projecting onto the instantaneous rest space of the 
fundamental observer, Ohm's law becomes J a = o~E a . 

We consider a model of the universe filled with a charged perfect fluid, which 
could be matter or radiation, characterised by the energy-momentum tensor Tp 3 = 
pu a u^ + ph alS ; the charge will give rise to an electromagnetic field, with T^f M = 
F a 1 F f3 ' y — jg a ^FjsF' yS . The energy-momentum conservation equation accounts for the 
interaction between the fluid and the field, so that [19 

Tf;p = F^Jfi ■ (9) 
Using the above definitions, the coupling term can be rewritten as 

F a/3 jp = qE a + F af3 J p = qE a + u a E p Jp + ff^JpR, . (10) 

We further assume the usual equation of state for a perfect fluid p = wp; then, the 
energy conservation equation for the fluid, m q T^ /3 ^ = u a F al3 jp, takes the form 

p + (l + w)Qp=-E l3 Jp, (11) 

and the momentum conservation equation h a pTp ' = h a pF^j^ takes the form 

(l + w)pu a + wD a p = qE a + e a ^J p B^ . (12) 

If the charged fluid is composed by non-relativistic matter, one defines q = een mat = 
—Pmat, where ee and n mat represent respectively the charge and the number density of 
the fluid particles, and m is their mass; if instead the radiation component is charged, 
one has q = e en ra( j oc p^J, with similar definitions. It appears from this picture that the 
introduction of a charged fluid breaks the isotropy of the universe, through the creation 
of currents and electromagnetic fields. Therefore, it is in principle possible to use the 
CMB as a measure of the isotropy of the universe to constrain the presence of an overall 
charge asymmetry. 

We must further consider the fact that the universe has a very high conductivity, 
during all the phases of its evolution. The conductivity varies in the different epochs 
of the evolution of the universe, depending on the number of charge carriers and on 
the dominant scattering process. Before e + e~ annihilation, the conductivity has been 
evaluated in [20], and is found to decay with temperature like o ~ T/(a 2 log(l/a)), 
where a is the fine structure constant. During the radiation dominated era, considering 
Thomson scattering and proceeding as in (21], one finds o = 2tt^^- ~ 3 x 10 13 sec -1 . 
During the matter dominated epoch, the dominant scattering process is again Thomson 
scattering of CMB photons with the residual free electrons, and one has o ~ 3 x 
10 10 sec _1 [22] . We perform our analysis in the ideal magnetohydrodynamics limit, 
for which the conductivity goes to infinity while the current remains finite ^H]- In 
the reference frame of comoving observers, corresponding to the reference frame of the 
fluid in this picture, Ohm's law takes the form J a = o~E a . Consequently, applying the 
magnetohydrodynamic limit, one finds that the electric field (but not the magnetic field) 
must go to zero |2S]- This ensures that the spatial current remains finite, avoiding the 
possibility of an instantaneous response of the fluid to the electromagnetic fields. Note 
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that this is valid because our analysis is performed in the reference frame of observers 
comoving with the fluid, characterised by the 4- velocity u a ; in a reference frame with 
4- velocity u a = u a — v a , the magnetohydrodynamic limit implies the presence of an 
electric field E a = —e^v^B 7 . 

Another way of rephrasing the magnetohydrodynamic limit is found by considering 
the time evolution of the electric and magnetic fields given by Maxwell's equations 
(0) and (|B|). The presence of a primordial, uniform magnetic field in a Friedmann- 
Robert son- Walker (FRW) universe is accounted for by assuming that the energy density 
of the field is a first order quantity in perturbation theory, and that the field itself is a 
small background component. We introduce the presence also of an electric field with 
the same characteristics. Then, the evolution of these fields at the lowest order is given 
by 

B a = - ^QB a , (13) 

E a = -J Q -^QE a . (14) 

The first of these equations suggests that the magnetic field varies on an approximative 
time-scale given by the Hubble time _1 . By substituting with Ohm's law J a = aE a 
in the second equation, one finds that the time-scale of variation of the electric field 
is instead of the order (a + (2/3)0) _1 . One further has that, for all the epochs of the 
evolution of the universe, cr/O ^> 1 (for example, o~/Q ~ 10 28 today). Therefore, while 
the magnetic field varies on a time-scale comparable to the Hubble time, the electric field 
varies on a much smaller time-scale, given by a -1 , and gets dissipated much quicker. 
The infinite conductivity limit provides an explanation of the reason why large scale 
magnetic fields, and not electric fields, are observed in the universe. 

In the case of a charged cosmic fluid which we are analysing, and in the reference 
frame of comoving observers, adopting the infinite conductivity limit reduces the form 
of the relevant equations to 

p+(l + w)0p = O, (15) 

(l + w)pii a + wD a p = e^J^B 1 , (16) 

for the energy and momentum conservation, and 

2u a B a =q (17) 

D a B a = (18) 

e^D^B 7 = J a - e al3 yB^ (19) 

hcpBP = (uj a0 + a a/3 -^eh a ^B p , (20) 

for Maxwell's equations. We have defined the charge density as q = een, both in the 
case of charged non-relativistic matter and radiation; therefore, using p mat cx n mat and 
p ra d oc n^, Eq. (fTKJ) implies the charge evolution q + Qq = 0. Equation (fTTjl shows 
that in a universe with infinite conductivity, the presence of a charge density implies 
the presence of a magnetic field, and induces vorticity in the metric. Therefore, the 
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analysis in a non-perturbative framework should be carried on in the context of tilted 
Bianchi universes, i.e. homogeneous and anisotropic universes in which the surfaces of 
homogeneity are not orthogonal to the matter flow. However, the goal of this analysis 
is to constrain the presence of charge in the universe by using the isotropy of the CMB. 
Consequently, the charge has to be considered a small perturbation, e — > 0, and we 
want to use linear perturbation theory on a FRW universe. The relevant equation at 
this purpose is (|17|): a constraint on the charge density can be derived from constraints 
on the magnetic field and the vorticity induced by it. We need to evaluate the time 
evolution of these two quantities. At the lowest order, the magnetic field evolution is 

B a + ^QB a = 0, (21) 

and the vorticity evolution is |2U 

u a + ~ecu a = -~e af37 D^. (22) 

Substituting with the momentum conservation equation (|16p. and using the identity 
e a01 D^D^f = -2fu a |21j, equation fl22D becomes 

u a + ( \ ~ w] Qu a = - 1 e a ^D^ Sfl J s B" . (23) 
\3 / 2p(l + w) 

The magnetic field evolution equation gives the usual scaling behaviour B a = 
Boa/a 2 , where B denotes the field amplitude today. The evolution for the vorticity 
is more involved. We know however that the charge density evolves like the number 
density of particles, q + Qq = 0: this equation can give us insight on the behaviour of 
vorticity. By deriving Eq. (|T7j) . and imposing the charge scaling, 2Co a B a -\-2ui a B a = — Qq, 
we find, from the magnetic field evolution 

B a u a = -~GB a u a . (24) 
o 

This identity is satisfied in the two cases (au a )' = 0, and (aw)' _L B. In the first case, 
the evolution of the vorticity is such that Co a + - ui a = 0; in the second case, one has 
that only the component parallel to the magnetic field satisfies io\\ + - u\\ = 0. In the 
following, we will always assume that uo a oc a -1 . Given that we want to constrain the 
charge using Eq. 1)17)1. we are in fact interested only in u\\. However, in the following we 
will set constraints on the vorticity vector independently, through the induced CMB 
anisotropies: our assumption, therefore, corresponds to postulating that the entire 
vorticity contribution to CMB anisotropies comes only from the u\\ component. A 
premise that is conservative, in the sense that it makes the final limit on the charge less 
tight. 

We now proceed to evaluate the bounds on the charge in two different 
configurations: uniformly and stochastically distributed charge. The strategy of limiting 
the charge asymmetry of the universe by using the observed isotropy of the FRW 
spacetime has already been adopted in two previous works: Orito and Yoshimura [2E] 
also used the measurement of CMB temperature fluctuations, instead Masso and Rota 
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[2*7] derived a bound using Nucleosynthesis. In both these references, however, it is 
claimed that the presence of a non-zero charge density would generate a large scale 
electric field in the universe: they modelled the universe as an insulating medium rather 
that a highly conducting one. Our work differs from both of them in this, we believe, 
crucial aspect. 

3. Uniform distribution of charge 

We first first focus on a uniform charge distribution in a homogeneous spacetime. The 
anisotropic expansion of the universe will leave an imprint on the propagation of light 
rays from the last-scattering surface until today. In 28 j tight constraints were derived 
on both uj and B for a general class of homogeneous, anisotropic models. It was found 
that 

B(t ) < 3 x l0- g tt 1/2 h Gauss , (25) 

u(t ) < 1(T 7 H , (26) 

where H = 100 /ikmsec _1 Mpc~ 1 is the Hubble constant today, and 1 — O is the 
curvature. The bound in Eq. ([SB*]) was found assuming the usual scaling for the vorticity 
in the matter era in the absence of sources: u a = uo a /d 2 , as can be derived from 
Eq. ("22]) • Therefore, we have to correct this bound accounting for our assumption for 
the evolution of vorticity sourced by the charge: u a = Uo a /a. This simply changes the 
limit by a factor of l/a(j] Tec ) ~ 10 3 : c^o < 10~ 4 H Q . With these constraints we obtain a 
limit on the charge asymmetry using the Schwarz inequality on Eq. ()17j) . 

q < 2\u(r )\\B(r )\ < 2.4 x 10~ 74 h 2 Q 1/2 GeV 3 . (27) 

It may be that the charge asymmetries corresponding to individual particle species 
are such as to cancel to some extent, yielding to a universe which is more neutral overall. 
However, should this not be the case, this last equation imply a maximal allowed charge 
Qx = q/n>x for each different constituent X, where n x is the number density of charged 
particles. For a flat universe (0 = 1), and assuming that the charge is due to a difference 
between the charges of the electron and the proton, we divide by the baryon density of 
the universe VLsh 2 = 0.02 to find 

g e _ p <10- 26 e. (28) 

If the charged particles are instead dominated by the dark matter constituents, 
we get gam ^ 4 x 10 _27 e m& m GeV -1 . For example, a conservative limit comes from 
neutralinos of mass of 1 TeV, gdm ^4x 10~ 24 e; considering instead a light dark matter 
particle with m dm ~ 10 MeV [21], one gets gdm ^ 4 x 10~ 30 e. If we assume that 
the charge imbalance is predominantly caused by photons, with n 7 = 421.84 cm -3 , we 
obtain g 7 < 10~ 35 e; for each species of neutrinos we obtain instead q u < 0.4 x 10 _34 e, 
from n u = 115.05 cm -3 . All these constraints are orders of magnitude more restrictive 
than the limits from laboratory experiments or astrophysical observations mentioned 
above, subject to the assumption that any charge asymmetries between species are not 
anti-correlated. 
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4. Stochastic distribution of charge 

One might expect that any mechanism of charge creation will be local, possibly a result 
of local violation of gauge invariance arising at some phase transition, or due to the 
escape of charge into extra dimensions. It would then be more appropriate to consider a 
stochastic distribution of charge asymmetry with, or without, a net imbalance in charge. 
The simplest assumption we can make, is that the magnetic field and the vorticity 
created by the charge density are two independent stochastic variables, with Gaussian 
distribution. Homogeneity and isotropy forces the first moment of the distributions to 
be zero and hence (q) = 0. We further assume that the magnetic field and vorticity 
induced have power law power spectra, B 2 (k) = Bk n and u 2 (k) = £1 k m , where the two 
spectral indices are different consequently to the independence assumption. We have 
then that: 



where we need to introduce an upper cutoff frequency k c in the magnetic field 
power spectrum, which accounts for the interaction of the magnetic field with the cosmic 
plasma at small scales; k c {r rcc ) ~ lOMpc -1 , from [3U|. For the following calculations, 
we do not need to assume a cutoff frequency for the vorticity. The factors (8ij — kikj) 
come from the divergence-free property of both the magnetic field and the vorticity [25J. 
We define B\ = (S*(x)5j(x)}|> and oj\ = (cj*(x)o; i (x))|A to be the energy densities of 
the magnetic field and the vorticity in a region of size A. 

We wish to estimate the mean fluctuation of charge in a region of size A: 





From Eq. (fTTj) one has 




(29) 



(30) 



< 9a ) - ^-t^^-st f(n, m, A, k c ) , 




(31) 



with 
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where T(a,x) denotes the incomplete gamma function (6.5.3 of |31|). 

Eq. (}3"T]) provides a limit for the charge asymmetry parameter q x = V\9a) /n>x 
today, as a function of the parameters characterising the magnetic field and the vorticity: 
amplitudes, spectral indices, cutoff frequency and coherence scale. 

The amplitude B\ of a cosmic magnetic field must satisfy stringent constraints, 
which have been found in previous works using completely independent methods. The 
most direct constraint on the strength of a magnetic field comes from the observation 
of Faraday rotation in radio sources: for example, in [32J it is shown that, for a cluster 
magnetic field, B\ < 10~ 6 Gauss on scales A ~ 1 kpc. If the charge imbalance, and 
consequently the magnetic field, was created before recombination, we can apply the 
limits which come from CMB observations (see for example [33]): i n [HI] limits where 
derived, by considering the anisotropies induced in the CMB by gravitational waves, 
created by the anisotropic stresses of a primordial magnetic field (see Eq. (28) of [34J, 
for A = 0.1/i -1 Mpc). Finally, in , a constraint was found on the amplitude of a 
magnetic field created before Nucleosynthesis, as a function of the spectral index n, by 
imposing the Nucleosynthesis bound on the energy density of the gravitational waves 
induced again by the field (Eq. (33) of |3S], also for A = 0.1/i _1 Mpc). 

We next need to find a constraint on u\, which is not present in the literature. 
Again, it can be found using the isotropy of the CMB. We need to evaluate the vector 
contribution to the CMB temperature anisotropy, and we use the formalism of [36J. 
Our vorticity variable Ui = ^e^io^ is related to the vector gauge invariant variable 
V c = where u' 1 ' is the vector part of the fluid velocity perturbation and B^ 

the vector metric perturbation, through relation uiy = aVcQuu]) where denotes 
the vector harmonic [2H]- Therefore, Ui = — B^)Qf\ For the purposes of this 

analysis, we neglect the finite thickness of the last scattering surface and work in the 
tight coupling limit, for which i^ 1 ) = v y = vb- In this limit, the dominant effect of 
vector perturbations is the creation of a dipole in the temperature anisotropy. From the 
general expression for the CMB power spectrum from vector perturbations given in [SB] 
we finally find 

0<= i^ r^f^o „><*), (3 2 ) 



7T a 2 (r rcc ) J (kr* 

where T* = tq — r rcc , and the factor a~ 2 (r rec ) accounts for the vorticity time evolution. 
This gives us 

£ _ . a ri 1 I 2 I 2 , ,2 * nm+l /oo\ 



2vr v r(^)r(^p) rec X T™+ 1 

The maximum value for the CMB temperature anisotropy is £ 2 C^/(27r) ~ 10~ 10 , and in 
order to constrain uj\ we fix the i dependence by considering the two extreme values 
i = 4, if m < -1, and t = 200, if m > -1. 

We can finally evaluate Eq. (|3*Tj). and derive an upper bound on the charge density 
in a region of size A, by applying the aforementioned constraints on B\ and the limit on 
oj\ just derived. The bounds are represented in fig. [TJ as a function of the spectral index 



Constraints on the Electrical Charge Asymmetry of the Universe 



10 




Figure 1. This figure shows the limits on the charge fraction parameter y which 
apply to a stochastic distribution of charge, as a function of the spectral index of the 
induced magnetic field, for A = 0.1/i _1 Mpc. The diagonally-shaded region represents 
the exclusion region derived from nucleosynthesis constraints |H5| : the vertically- 
shaded region is excluded by CMB anisotropies [HI]- Magnetic field constraints from 
clusters \->2\ give a less restrictive bound, which is not shown in the figure. If the 
process generating the charge is causal, then n > 2. 



n of the magnetic field. For every value of n, we have maximised (|31|) with respect to 
the vorticity spectral index m, on which we have no theoretical prediction. Note that 
the quantity shown in the y-axis, 

(34) 

e \U B m x J 

is the charge fraction parameter 'normalised' to the case of 'baryonic' charge asymmetric 
matter (as in Eq. ([28))) . If the process generating the charge and the magnetic field is 
causal, e.g. it takes place during a phase of standard Friedmann expansion, then the 
spectral index must satisfy n > 2 [3*7j . In this case, the limits on qx are more stringent 
than in the case in which the generation occurs during inflation. 



5. Discussion 



In this paper we have derived a cosmological constraint on the presence of a non-zero 
electric charge in the universe. A conservation law associated with a long range force, 
such as charge conservation with the electromagnetic force, is usually believed to be 
exact as a result of gauge invariance. Conservation of charge appears to hold in all 
particle decays, and there are strong experimental constraints on the charge of particles 
which are predicted to be neutral by the standard model. However, even though electric 
charge conservation is well established on Earth, this may not imply directly the overall 
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neutrality of the universe: to be able to draw this conclusion, one would have to assume 
in addition that charge has to be conserved on all scales and during the entire evolution 
of the universe. This may not be the case and the determination of a cosmological 
constraint on the charge asymmetry of the universe is of conceptual importance. 

There are a host of proposals that lead to the generation of a charge asymmetry in 
the universe. In the context of brane world scenarios, charge non-conservation arises in 
the form of charge leakage into extra dimensions; in varying speed of light theories, it 
can be a consequence of the variation of the fine structure constant; allowing for some 
modification of the standard model, it can arise if the gauge symmetry is temporarily 
broken during a phase transition taking place in the early universe, or if one imposes 
a small, non-zero mass to the photon. In our analysis, we did not consider a complete 
model in which the charge asymmetry originates: by doing so, our constraints gain in 
generality. The assumption we made, that charge is conserved in the period in which our 
constraint is established, points toward a model in which the charge is created during 
a primordial phase transition leading to a transitory breaking of the electromagnetic 
gauge symmetry. 

We have obtained our constraints on the overall charge imbalance of the universe 
in two different cases: a uniform distribution of charge, and a distribution characterised 
by stochastic fluctuations. In order to derive our results, we made the assumptions that 
the universe is a good conductor, and that the charge is a first order perturbation 
in a background FRW model; moreover, we generalised the scaling of the charge- 
induced vorticity as the inverse of the scale factor. These are the only assumptions 
which affect the uniform distribution limit; in the case of the stochastic distribution, 
we further assumed that vorticity and magnetic field are two independent, Gaussian 
stochastic variables. A possible extension of our analysis would be to go beyond the 
linear calculation, and evaluate the full effect that the presence of a non-zero charge has 
on the background dynamics of the spacetime. 

In the case of a uniform distribution of charge, it is remarkable that our cosmological 
limit, once translated in terms of constraints for the single charge carriers, gives results 
which are orders of magnitude stronger than the ones derived by terrestrial experiments 
or astrophysical observations. If the charge is distributed stochastically, the limits are 
less stringent: this is a consequence of the fact that the CMB bound on stochastic 
vorticity is less tight than in the uniform case, and in addition we have maximised 
it with respect to the vorticity spectral index m. However, we still find interesting 
constraints for high values of the magnetic field spectral index: these would apply, in 
the case of a causally created charge asymmetry, for example during a phase transition. 
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